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Abstract 



I review the notion of the quark-hadron duality from the modern perspective. 
Both, the theoretical foundation and practical applications are discussed. The 
proper theoretical framework in which the problem can be formulated and treated 
is Wilson's operator product expansion (OPE). Two models developed for the de- 
scription of duality violations are considered in some detail: one is instanton-based, 
^ another resonance-based. The mechanisms they represent are complementary. Al- 

though both models are rather primitive (their largest virtue is their simplicity) they 
Ph. hopefully capture important features of the phenomenon. Being open for improve- 

ments, they can be used "as is" for orientation in the studies of duality violations 



in the processes of practical interest. 
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27-29 March 2000. To be published in the Boris Ioffe Festschrift At the Frontier of Particle 
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1 Introduction 

Quantum chromodynamics (QCD) is a very strange theory. All theoretical calcula- 
tions are done in terms of quarks and gluons. At the same time, quarks and gluons 
are never detected experimentally. What is actually produced and detected in exper- 
imental devices are hadrons: pions, kaons, protons, etc. The quark-hadron duality 
allows one, under certain circumstances, to bridge the gap between the theoretical 
predictions and experimentally observable quantities. The idea was first formulated 
at the dawn of the QCD era by Poggio, Quinn and Weinberg |I|], who suggested 
that certain inclusive hadronic cross sections at high energies, being appropriately 
averaged over an energy range, had to (approximately) coincide with the cross sec- 
tions one could calculate in the quark-gluon perturbation theory. The uses of this 
theoretical construction are countless: the e + e~ annihilation, deep inelastic scatter- 
ing, hadronic r decays, inclusive decays of heavy quarks, physics at the Z peak, to 
name just a few. In spite of the dramatic developments in QCD in the subsequent 
two decades the notion of the quark-hadron duality remained vague, essentially at 
the 1976 level, with the very basic questions unanswered. These questions are: 



• What energy is considered to be high enough for the quark-hadron duality to 
set in, and what accuracy is to be expected? 

• What weight function is appropriate for the averaging of the experimental cross 
sections? 

• If the theoretical prediction includes only perturbation theory, should one limit 
oneself to some particular order in the the a s series? 

• Do we have to include known nonperturbative effects (e.g. condensates) in the 
theoretical prediction? 

• Given a definition of the quark-hadron duality, can one estimate deviations 
from duality and how? 



Of course, not all of these questions are independent. For instance, answering 
the last question, one will simultaneously learn the boundary energy. 

Systematic explorations of these and related issues started in earnest in 1994 [Q . 
As in many other instances, this was dictated by practical needs. Previously, the 
accuracy of the experimental data on hard inclusive processes was rather modest, so 
that the Poggio-Quinn- Weinberg prescription was good enough. By 1994 the data, 
mostly associated with the b quark physics, became so precise and the questions 
raised so acute, that a much better theoretical understanding became imperative. 
Probably, the most clear-cut example is the problem of the B semileptonic branching 
ratio ||: theoretical expectations obtained in the quark-gluon language exceed the 
measured number by 10 to 20%. Possible deviation from duality is suspected to be 
a major source of theoretical uncertainties. If one could reliably rule out duality 



violations at this level, the conclusion of the interference of new physics would ensue 
(provided the experimental data stay intact, of course). The stakes are quite high. 

It is fair to say that (short of the full solution of QCD) understanding and 
controlling the accuracy of the quark-hadron duality is one of the most important 
and challenging problems for the QCD practitioners today. In this issue one cannot 
expect help from lattices. The duality violation is a phenomenon inseparable from 
the Minkowskian kinematics; numerical Euclidean approaches, such as lattice QCD, 
have nothing to say on this issue. Analytic methods are needed. 

In this review I will summarize the results of an investigation |l [10|, which spans 



over six years and is not yet complete. I will discuss the modern formulation of the 
problem based on Wilson's operator product expansion (OPE). This formulation is 
solid and unambiguous. Then I will review models which were designed to give us 
a certain idea of (and a degree of control over) the deviations from duality. There 
are two classes of such models: instanton-based and resonance-based. They are 
complementary to each other; I will discuss both classes. Finally, I will present 
sample applications in the processes of the current interest, such as the hadronic r 
decays. 

2 The Quark-Hadron Duality: What Does That 
Mean? 

Let us consider an idealized theory - QCD, with two massless quarks, u and d. 
We are interested in the total hadronic cross section of the e + e~ annihilation. The 
"photon" in our theory is idealized too, its coupling to the quark current has the 

form 

u-j^u - d-j^d 

(In fact, this is the isovector part of the actual electromagnetic current). We define 
the two-point function LV 

IV = t J e lllx d A x(0\T{Ux)Ji(0)}\0) . (2) 

Here q is the total momentum of the quark-antiquark pair. Due to the current 
conservation LV is transversal, 

IV = {q^q v - q 2 g flu )Tl(q 2 ) . (3) 

The experimentally observable quantity is the imaginary part of II(g 2 ) at positive 
values of q 2 (i.e. above the physical threshold of the hadron production), the spectral 
density, 

127T 

p(s) = —ImU(s) , s = q 2 . (4) 




Figure 1: The one- loop graph determining the polarization operator and the spectral 
density in the leading (parton) approximation. The "photon" momentum is denoted 
by q. 

Up to a normalization, the expression above coincides with the cross section of e + e~ 
annihilation into hadrons measured in the units a(e + e~ — > fi + [i~), the famous ratio 
R. 

Theoretically one can calculate II (q 2 ) in the deep Euclidean domain, at negative 
q 2 . For instance, from the free quark loop of Fig. 1 one gets 

n(g 2 )-- I ^lnQ 2 Q 2 = -q 2 . (5) 

Performing the analytic continuation to the Minkowski domain and taking the imag- 
inary part one arrives at 

P0)thcor -> 1 , S -> OO . (6) 

This is the spectral density in the theory with free quarks, i.e. a s = 0. There 
are various corrections to the free quark result (|^). The perturbative gluon ex- 
changes give rise to the a s (Q 2 ) series. Nonperturbative (power) corrections come 
from the quark and gluon condensates and from other sources, e.g. the small-size 
instantons. A systematic method of handling the theoretical calculations of H{Q 2 ) 
in the deep Euclidean domain is provided by Wilson's operator product expansion 
(OPE) |l IJ. Essentially, this is a bookkeeping procedure: one consistently separates 
the short-distance contributions (i.e. those coming from distances < /i -1 ) from the 
large distance contributions (i.e. those coming from distances > /i -1 ). Here ji is 
a theoretical parameter (usually referred to as the normalization point) separating 
the two domains. The choice of /i is a matter of convenience - observable quantities 
do not depend on it. 

The short-distance contributions determine the coefficients C n (q) of OPE , 

D{q 2 ) =Y,C n {q^)(OM) , D(Q 2 ) = -(An 2 )Q 2 (dU/dQ 2 ) . (7) 

n 

The normalization point \x is indicated explicitly. The sum in Eq. (|7|) runs over all 
possible Lorentz and gauge invariant local operators built from the gluon and quark 
fields. The operator of the lowest (zero) dimension is the unit operator I, followed 
by the gluon condensate G 2 , of dimension four. The four-quark condensate gives 
an example of dimension-six operators. 

At short distances QCD is well described by the quark-gluon perturbation theory. 
Therefore, as a first approximation, it is reasonable to calculate C n perturbatively, 



in the form of expansion in a s (Q) ~ (InQ) -1 . This certainly does not mean that 
the coefficients C n are free from nonperturbative (nonlogarithmic) terms. The latter 
may and do appear in C n 's; they are of the type ~ Q~~' where 7 is a positive number, 
not necessarily an integer. Such terms in C n 's are generated, for instance, by the 
small-size instantons. Another source of the power terms in C n , of a technical rather 
than dynamical nature, is the normalization point p: in calculating C n (fi) one must 
remove all soft contributions with off-shellness less than p. 

The condensate terms in Eq. (|7|) give rise to corrections of the type (Aq CD /Q) n 
where n is an integer > 4, the (normal) dimension of the operator O n , modulo 
logarithms associated with the anomalous dimensions. Aqcd is the scale parameter 
of QCD (sometimes, I will drop the subscript "QCD" for brevity) entering through 
the matrix elements of CVs. 

If one could calculate II (Q 2 ) in the Euclidean domain exactly, one could analyt- 
ically continue the result to the Minkowski domain, and then take the imaginary 
part. The spectral density p(s) t heor obtained in this way would present the exact 
theoretical prediction for the measurable hadronic cross section. There would be no 
need for duality. 

In practice, our calculation of IT(Q 2 ) is approximate, for many reasons. First, 
nobody is able to calculate the infinite a s (Q 2 ) series for the coefficient functions, let 
alone the infinite condensate series. Both have to be truncated at some finite order. 
A few lowest-dimension condensates that can be captured, are known approximately. 
The best we can do is analytically continue the truncated theoretical expression, term 
by term, from positive to negative Q 2 . For each term in the expansion the imaginary 
part at positive q 2 (negative Q 2 ) is well-defined. We assemble them together and 
declare the corresponding p(s) t hcor to be dual to the hadronic cross section p(s) exp . 
In the given context "dual" means equal. 

Let me elucidate this point in more detail. Assume that Tl(Q 2 ) is calculated 
through a 2 and 1/Q A , while the terms a 3 and 1/Q 6 (with possible logarithms) are 
dropped. Then the theoretical quark-gluon spectral density, obtained as described 
above, is expected to coincide with p(s) cxp , with the uncertainty of order 0[(a s (s)) 3 ] 
and 0(1/ s 3 ). The uncertainty in the theoretical prediction of this order of magnitude 
is natural since terms of this order are neglected in the theoretical calculation of 
H(Q 2 ). If the coincidence in this corridor does take place, we say that the quark- 
gluon prediction is dual to the hadronic spectral density. If there are deviations 
going beyond the natural uncertainty, we call them violations of duality. Needless to 
say that, once our calculation of Tl(Q 2 ) becomes more precise, the definition of the 
"natural uncertainty" in p(s)theor changes accordingly. 

This is the most clear-cut definition I can suggest. From the formal standpoint, 
it connects the duality violation issue with that of analytic continuation from the 
Euclidean to Minkowski domain. Negligibly small corrections (legitimately) omitted 
in the Euclidean calculations may and do get enhanced in Minkowski. 

Before I proceed to explain, from the physical standpoint, where the violations 
of duality come from, and their pattern, I would like to make a remark important in 



the conceptual aspect. The necessity of truncation of the a s and condensate series 
is not just due to our practical inability to calculate high-order terms. (For the 
vast majority of the theorists, myself including, "high-order" begins at the next- 
to-next-to-leading level.) Assume we have a "sorcerer's stone" which would allow 
us to exactly calculate any term in the expansion we want. Still, we would not be 
able to find n((5 2 ) CX act because the both series are factorially divergent. The a s 
series has at least two known sources of the factorial behavior: first, the number of 



the Feynman graphs grows factorially at high orders jE^; second, there are graphs 
with renormalons, (of which we will not bother in what follows since the infrared 
renormalons - the only ones which are potentially dangerous - are totally eliminated 
by the introduction of the normalization point \x ||13|| ). The condensate series is 
divergent too. The factorial divergence of the condensate series is studied even to 
a lesser extent than that of the a s series. The only fact considered to be firmly 
established [0, |J is the divergence per se. It is not Borel-summable, which sets a 
limit of the theoretical accuracy. Including more and more terms in the series would 
not help, and even an optimal truncation would leave a gap between n(Q 2 ) t heor and 

II(Q ) exact- 

3 Where the Duality Violations Come From? 

Theoretical calculations of inclusive processes in QCD are performed through Wil- 
son's OPE in the Euclidean domain. Therefore, to understand what is included in 
such calculations and what is left out, one should have a clear picture of limitations 
of OPE. 

Let us return to the consideration of the correlation function of two currents, 

(0|T{J,(x)jt(0)}|0), 

which determines the polarization operator in Eq. (|2]), at negative (Euclidean) x 2 . 
Wilson's OPE is nothing but an expansion of (0|T{ J M (x) Jt(0) }|0) in singularities at 
x 2 = 0. It properly captures all terms of the type 1/x 6 times logarithms, or lnx 2 , or 
x 2 \nx 2 , and so on. Every term of this type is represented in the condensate series 
provided that the calculation is carried out to a sufficient order. 

Let me note in passing that limiting oneself to singular terms (and, in particular, 
to the leading singular term) is a reasonable approximation in the theories in which 
the asymptotic conformal regime sets in at short distances in the power-like manner. 
In fact, OPE in its original formulation was designed by Wilson for applications 
in such theories (see the first work in Ref. |11|). In QCD the approach to the 
asymptotically free regime is very slow - logarithmic. That's why the precision of 
the leading (parton) approximation is sufficient only for very rough estimates, and 
that's why the high-order terms, both logarithmic and power, must be kept. The 
time when the QCD practitioners would be satisfied by the leading, or even the first 
subleading approximation, is long gone. 



It is clear that the function (0|T{J M (x) Jt(0)}|0) is not fully determined by its 
singularities at x 2 = 0. Generally speaking, one may have, additionally, isolated 
singularities at finite x 2 , or a singularity at infinity, which are not reflected in the 
truncated Wilson's OPE. Consider, say, a singularity at finite x 2 of the form 

(8) 



x 2 + p 



The expansion of this function (truncated at any order) generates derivatives of 
S(q 2 ) in the momentum space; therefore, it has no impact on the OPE expansion 
at large Q 2 . This contribution is clearly missing. The Fourier transform of Eq. ([8]) 
at large Euclidean Q 2 falls off as exp(— Qp). Thus, this term is smaller than any 
of the terms in the condensate expansion. We should not forget, however, that our 
final goal is a prediction in the Minkowski domain. Upon analytic continuation, the 
exponentially small term exp(— Qp) looses its suppression and becomes oscillating, 
sin(— Ep), where E stands for the total energy, E = yfq 1 . If it were not for the 
power suppression in the pre-exponent, such exponential/oscillating terms would be 
a total disaster. Since they are not seen in OPE, any prediction for the inclusive cross 
sections made through OPE (which is equivalent, in practice, to perturbation theory 
plus a few condensates) would be grossly wrong. There would be no quark-hadron 
duality at all. All calculations of the hard processes - from the total hadronic cross 
sections in e + e~ annihilation, to jet physics, to the heavy quark decays - would be 
valid roughly up to a factor of two, no matter how large the energy release is. 

Fortunately, one avoids the disaster due to the fact that the duality-violating 
exponential/oscillating terms are, in fact, suppressed as E~ K sin(-Ep) where k is 
a positive index, which depends on the process under consideration and is typically 
rather large. This justifies theoretical predictions based on a few first terms in 
the condensate expansion. Needless to say that determining k is one of the most 
important tasks in the issue of the quark-hadron duality. 

Singularities at x 2 — > oo also lead to the exponential/oscillating terms of a 
somewhat different form, of the type exp(— Q 2 p 2 ) in Euclidean. In the Minkowski 
domain one gets an oscillating function of the type E~ v sin(—E 2 p 2 ), where 77 is 
another positive index. In one and the same process, one can expect both duality- 
violating components to show up. Generally speaking, the indices k and rj are 
unrelated; at least, at the moment we do not see any obvious relation between 
them. (The coincidence of k and rj in the practically important problem of the 
inclusive hadronic r decays, k = T] = 6, seems to be accidental, see Sec. 8). 

How can one isolate the singularities at x 2 = from others diagrammatically? 
To answer this question it is convenient to pass to the momentum space. In the 
leading approximation the polarization operator H(Q 2 ) is presented by the graph of 
Fig. 1. The large Euclidean momentum q (remember, Q 2 — > 00) flows in the photon 
line on the left, propagates through the fermion lines, and leaves the graph through 
the photon line on the right. The virtual momenta of the fermion lines typically 
scale as q. The first a s correction is presented in Fig. 2. In this graph the virtual 
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Figure 2: The one-gluon correction in the the polarization operator. The gluon 
momentum is denoted by k. 

momenta of all lines in the loops, the fermion and the gluon, scale as q. Thus, the 
diagram of Fig. 2 determines the leading logarithmic correction to the parton result 
(|)or(§. 

Although the leading contribution to the integral comes from the domain where 
all virtual momenta are proportional to q, there are nonvanishing contributions 
from other kinematical domains. For instance, one can consider a "corner" where 
the gluon virtual momenta k in Fig. 2 is small, \k\ < fi, and does not scale with q. 
Certainly, at small k the gluon propagator is not given by perturbation theory. The 
gluon line must be cut. This corresponds to the gluon condensate term |14j in the 
condensate expansion. 

In general, any term in the condensate expansion can be interpreted in this way, 
in terms of factorization in the momentum space. The large external momentum q is 
transmitted through one or several hard lines - their momenta scale with q. This is 
the definition of "hardness." The remainder of the graph is a soft part, which factors 
out and gives rise to gluon, quark and mixed condensates. The virtual momenta in 
this part of the graph do not scale with q, they are limited by a fixed parameter \i 
(this is the definition of "softness"). 

The hard part of the graph is responsible for the coefficient functions. The fact 
that not all lines in the given graph are hard results in the power suppression of 
the corresponding coefficient function. Letting more and more lines to be soft, one 
obtains consecutive terms in the condensate expansion. What is missing? 

It is conceivable that the number of lines through which q is transmitted becomes 
so large, that though Q — » oo, neither of the lines is hard. Of course, in this case 
one cannot speak of the separate lines in the Feynman graphs. It would be more 
relevant to say that the external momentum q is transmitted from the incoming to 
the outgoing photon through a coherent soft field fluctuation, see Fig. 3. An example 
is provided, for instance, by a fixed-size instanton. It is clear that this mechanism 
is conceptually related to the truncated tail of the condensate series. Indeed, as 
one proceeds to higher condensates, more lines become soft. Eventually we arrive 
at the situation when all lines are soft. Mathematically, the exponential/oscillating 
contribution is related to the factorial divergence of the condensate series |4]]. This 
is the usual story. The exponentially small terms in Euclidean convert into an 
oscillating function in Minkowski. 

I will use the fixed-size instantons for the purpose of modeling this mechanism 
of duality violations. The observation that "soft" instantons generate an oscillating 





Figure 3: Transmitting a large external momentum through a soft field. 

component ascends to Refs. fl5| , |T6|j . By no means I imply that the instantons 
are the dominant soft fields in the QCD vacuum. True, there are models in which 
they are assumed to be dominant, the so-called instanton liquid models |T7|]. Their 
status in the range of questions I am interested in (the duality violations) has yet 
to be clarified. I will use instantons for the purpose of orientation, in hope that at 
least some features of the results obtained in this way will be more general than the 
model itself. 



4 Model or Theory? 

The topic I address - the quark-hadron duality violations - has a unique status. 
By definition, one cannot build an exhaustive theory of the duality violations based 
on the first principles. Indeed, assuming there is a certain dynamical mechanism 
(which goes beyond perturbation theory and condensates) for which such a theory 
exists, one will immediately include the corresponding component in the theoretical 
calculation. The reference quantity, n(Q 2 ) t h e0 r, will be redefined accordingly. After 
the analytic continuation to Minkowski, this will lead, in turn, to a new theoretical 
spectral density to be used as a reference p(s)theor in the duality relation. 

Thus, by the very nature of the problem, it is bound to be treated in models 
of various degrees of fundamentality and reliability. This is because the duality 
violation parametrizes our ignorance. Ideally, the models one should aim at must 
have a clear physical interpretation, and must be tested, in their key features, against 
experimental data. This will guarantee a certain degree of confidence when these 
models are applied to the estimates of the duality violations in the processes and 
kinematical conditions where they had not been tested. 



5 The Physical Picture Behind the Duality 

Before delving into technical details I will describe the phenomenon from a slightly 
different perspective. The quark-hadron duality takes place in those processes where 
one can isolate two stages in the process under consideration, occurring at two 
distinct scales. A basic transition involving quarks (gluons) must typically occur at 
a short scale regulated by external parameters such as Q, tuq, etc. For instance, 
in the e + e~ annihilation the basic transition is the conversion of the virtual 7 into 
uu or dd. Then, at the second stage, the quarks (gluons) materialize in the form 
of hadrons, at a much larger scale. In the appropriate frame, the first time scale 
is of order 1/Q while the second of order Q/A 2 . By that time, the original quarks 
are far away from each other - a residual interaction cannot significantly alter the 
transition cross section which was "decided" at the first (quark-gluon) stage. 

The duality violations are due to (i) rare atypical events, when the basic quark 
transition occurs at large rather than short distances; (ii) residual interactions oc- 
curring at large distances between the quarks produced at short distances. In the 
first case appropriate (Euclidean) correlation functions develop singularities at finite 
x 2 , while the second mechanism is correlated with the x 2 — ► 00 behavior. 

In both cases the duality violating component follows the pattern I have dis- 
cussed above - exponential in Euclidean and oscillating in Minkowski. Three distinct 
regimes were identified and considered in the literature so far: 

• (i) Finite-distance singularities 

s- K / 2 sin(v^) ; (9) 

• (ii) Infinite-distance singularities (N c = 00) 

s- v/2 sin(s) ; (10) 

• (iii) Infinite-distance singularities (N c large but finite, s — >■ 00) 

exp (—as) sin(s), a = O I — J <1. (11) 

These regimes are not mutually exclusive - in concrete processes one may expect 
the duality violating component to be a combination of (i) and (ii), or (i) and (iii). 
From the theoretical standpoint it is quite difficult to consistently define the duality 
violating component of the type (3). An operational definition I might suggest is as 
follows: Start from the limit N c = 00 and identify the component of the type (2). 
Follow its evolution as N c becomes large but finite. 

6 An Instanton-Based Model 

The basic features of the formalism [[| to be used below to model exponential/oscil- 
lating terms (which present deviations from duality) are as follows: 
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(i) One considers quarks propagating in the instanton background field. The 
instanton size p is assumed to be fixed. Alternatively, one can say that an effective 
instanton measure has a ^-function peak in p. 

(ii) For the light quarks, most instanton amplitudes are suppressed by powers 
of the light quark masses. This suppression is due to the fermion zero modes of 
the Dirac operator, occurring in the instanton-like backgrounds. We ignore these 
factors, as well as all other pre-exponential overall factors coming from the instanton 
measure. We will only trace dependences on large momenta relevant to the problems 
under consideration (Q in e + e~ annihilation, the heavy quark mass ttiq in the case 
of the heavy quark decays, and so on). 

(iii) We ignore all singularities of the correlation functions at x = 0. The cor- 
responding contributions are associated with the power terms in the condensate 
expansions. The instantons-based models are presumably not precise enough to 
properly capture the condensates. We will isolate and calculate only those contribu- 
tions that come from the finite distance singularities at x 2 = — p 2 . In the momentum 
space they produce the exponential/oscillating terms sought for. 

It is seen that our evaluation of deviations from duality is based on the most 
general aspects of the instanton formalism, and, in essence, does not depend on 
details. 

The advantage of the model is its simplicity. The only formula we will need is 

A - , 1 R *• = * m- 2 **-.W») (12) 

(x 2 + P 2 y i» v 2 J p 2 "- 4 v ; 

in the Euclidean domain, which implies that in the Minkowski domain, at large q 2 , 

/i 
d A x—^ — e iqx oc s^"4 sin( v / sp - 5) , s = q 2 . (13) 

V ' P ) 

Here K is the McDonald function, and S is a constant phase which is of no concern 
to us here. 

To explain how it works it seems best to consider concrete examples. Let us 
start from e + e~ annihilation. The polarization operator defined in Eq. (|2|) is given 
by the graph of Fig. 1. In Sec. 2 we evaluated this graph for free quarks and found 
that p(s) tends to a constant at asymptotically large energies. Now, we evaluate 
the very same graph using the quark Green functions in the background instanton 
field, rather than the free quark Green functions. For massless quarks the Green 
functions G; ns t are known exactly |lj| , but we will not need the full expression. We 



will only need to know that Ginst is a sum of terms of the following structure 

Ginst ^' y) = [(x-y)2]2 [(a;-z)2 + p2 ] £ l[(?/ _ z) 2 +p 2 ] £ 2 G > ( 14 ) 

where z is the instanton center, and 

ti? = l or ^, h + i 2 = 2, (15) 
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and the numerator G is a polynomial oix,y. As was explained above, the singularity 
of G inst (x,y) at x — y = is irrelevant for the exponential/oscillating terms. Of 
importance are the singularities in the complex plane coming from the second factor 



in Eq. (14). The integrals which one has to take can (and must) be evaluated at 
a saddle point; a simple analysis || shows that at the saddle point the instanton 
center z is exactly in the middle between x and y. Then the relevant singularities 
of the quark correlation functions are at (x — y) 2 = — 4p 2 (see Eq. ([14])). At the 
singularity the first factor l/(x — y) 4 (which is singular at the origin) can be replaced 
by l/(16p 4 ) and then safely omitted together with all other prefactors. 

We can proceed now to the calculation of the exponential/oscillating component 
of the polarization operator IT^. The polarization operator is the product of two 
Green functions (0); therefore, at large Euclidean momenta 



1 
'z-, 



V K J dX€ d ^[ {x ^ z y +p 2 ][{y _ z) 2 +p 2 ] 3 

oc f d 4 xe iqx — |— . x I d 4 ze^ z . - - -. 
J x 1 + p l J (z 2 + p 2 ) 6 

ex K 1 (Qp)K_ 1 (Qp)^^ex V (-2Qp). (16) 

Note that once the integration over the instanton center is carried out, the integral 
factorizes. 

Equation fliTf ) implies that the exponential component of n(g 2 ), defined in Eq. 
(i),is 

An oc ij exp(-2Qp) , (17) 

which implies, in turn, that at large E = y/s the oscillating component of the 
spectral density is 

Ap(s)<x-^sm(2Ep). (18) 



Equation ( |T8"D reproduces the high-energy asymptotics of the exact result [IS] for 
the polarization operator in the one-instanton approximation. 

Does the s _3//2 fall off of the oscillating (duality violating) component make 
sense? I will confront theoretical expectations with experiment in Sec. 11. Here I 
just refer the reader to Fig. 11 (see the dashed curve), postponing a more detailed 
discussion till after both, the instanton-based and the resonance-based models, are 
considered. 

The next example to be analyzed is the total hadronic r width. This exercise is 
quite similar to previous exercises. The corresponding transition operator is depicted 
in Fig. 4 (its imaginary part at p 2 = M 2 is proportional to the width of the hadronic 
r decay). The neutrino Green function clearly does not "feel" the background gluon 
field, so that the neutrino propagator is that of a free fermion. The same is valid 
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Figure 4: The transition operator T determining the hadronic r width, Thadrty) = 

M- 1 Im(r|f|r). 

for the r lines. Therefore, we immediately conclude that 

Ar(r -» 1/ + hadr.) oc — sin(2M r p) , (19) 

cf. Eq. (pf). The asymptotic (parton- model) prediction in this case is T(r — > 
z/ + hadr.) oc M%, so that the oscillating component in R T scales as 

Ai2 T oc-^sin(2M r p), (20) 

where 

__ r(r~ -► i/ T + hadrons) 

^r = TTT-I — TT • (21) 

Note that the pre-exponential suppression factor in Ai? r is significantly stronger 
than in Ai?(e + e~), namely, M~ 6 vs. E~ 3 . Of course, to make quantitative state- 
ments it is not sufficient to establish the scaling laws; one needs absolute normal- 
izations. Here, we are basically in uncharted waters. At best, we have ||] some 
educated guesses, which, if true, imply that AR T /R T ^ 5%. This estimate is rather 
close to what one obtains for duality violations in r in the resonance-based model, 
see Sec. 8. 

In Sec. 1 I mentioned that the current interest to the problem of the quark- 
hadron duality was driven, to a large extent, by a significant progress in the experi- 
ments on the inclusive heavy flavor decays. Theoretically, they can be treated along 
the same lines as e + e~ annihilation or r decays. The distinctions are technical. Let 
us start from the total semi-leptonic width of the b flavored hadrons. At the quark 
level the process is described by the transition 



where q is either u or c quark, and £ stands for the electron, muon or r lepton. We 
will first neglect the masses of the final fermions; this is an excellent approximation 
for q = u and C. = e, \x. The impact of the final quark (lepton) masses will be 
considered later. 

The relevant transition operator is depicted in Fig. 5. The differences compared 
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Figure 5: The transition operator T determining the total semileptonic width of B 
mesons. 

to the case of the r decay are as follows: 

(i) There is one (rather than two) light-quark line carrying a large external 
momentum; 

(ii) Unlike the r lines, the h quark lines do experience interactions with the 
background gluon field. 

The analysis of the exponential/oscillating component in the heavy flavor decays 
combines standard elements of the heavy quark expansion (for a review see e.g. 
|p0|), and the instanton calculus (a pedagogical review can be found in |2l]]). It is 
convenient to choose the rest frame of the heavy hadron at hand, and single out the 
large "mechanical" part in the x dependence of the heavy quark field, 

Q{x) = e- im ^Q(x) . 

Then the total width is proportional to the imaginary part of the transition operator, 

T = 1 ±-(H Q \f\H Q ), (22) 

M Hq 

where 

f = i I Q(x)G(x, y)Q(y)D(x - y)e im ^ x °- yo) d 4 (x - y)d 4 z , (23) 



G(x,y) is the light quark Green function, D(x — y) describes the propagation of 
colorless objects (the lepton pair in the case at hand), while z is the instanton 
center. The subscript marks the time component. 

Both, the light quark Green function G(x, y) and the heavy quark fields Q are to 
be considered in the background gluon field. Taken separately, they are not gauge 
invariant; the product in Eq. fl23|) is. In the leading order in the heavy quark 
expansion, the heavy quarks propagate only in time; therefore, 

Q(x , x) = Te'fo Ao{T ' g)dT Q(0, x) = U(x)Q(0, x) . (24) 

It is convenient (although not necessary) to impose the condition that at large dis- 
tances from the instanton center the quark propagation becomes free. This condition 
implies that the singular gauge is used. An explicit expression for U(x) can be found 
in this gauge; it is rather cumbersome, and I will not quote it here, referring the 
reader to the original publication || . At the saddle point (which again corresponds 
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to the instanton situated exactly in the middle between the points, x and y, i.e. 
z = (1/2) (a; + y)) the product U~ l (x)...U{y) in Eq. (|2~3|) reduces to unity, at least 
in the part which is singular at (x — y) 2 = —4p 2 . Note, that the matrices converting 
the nonsingular-gauge Green function G(x, y) to the singular gauge can be ignored 
too. This implies that one may continue to use Eq. flT3|). 

Thus, the heavy quarks decouple from the instanton background in the calcula- 
tion of the exponential/oscillating terms. This is the consequence of the fact that we 
exploit the heavy quark expansion and limit ourselves to the leading in 1/rriQ terms. 
In the subleading terms this decoupling does not necessarily take place. Generally 
speaking, the replacement of the z integral by the value of the integrand at the 
saddle point is not warranted in the next-to-leading orders. This effect will not be 
further pursued, however. 

Summarizing, the heavy quark expansion makes the heavy quarks sterile with 
respect to the duality violating component. This is certainly not counterintuitive. 
The duality violating component in the heavy quark inclusive decays originates from 
the light-quark propagators. Its calculation reduces to Eq. fllil) and the routine 
outlined after this equation (plus the basic formula (0)). Starting from Eq. ( |23|) 
we arrive at 

f oc \e~ 2m Q p (25) 

in the Euclidean domain, which results in 

Ar s i oc -^r sm(2m Q p) . (26) 

Since the total width scales as T oc ttiq, the oscillating component of the semileptonic 
branching ratio is obviously suppressed by tuq, 

ABt(Hq — > iv + hadrons) oc — 5- sin(2m Q p) . (27) 

It is easy to generalize this formula to include an arbitrary number of the light 
quarks in the final state. Each extra light quark adds two powers of rnq in the 
numerator. Thus, in the total nonleptonic branching ratio, with three light quarks 
in the final states, 

ABt(Hq — > light hadrons) oc — 7- sin(2mgp) . (28) 

m Q 

For the sake of completeness I will also mention the radiative decays of the type 
6 ^ s + 7. Assuming that these decays are induced by local operators^ 

fov(l + 7 5 )sF jUi ,, 



2 In actuality this is true only for a part of the amplitude. 
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where F^ u is the photon field strength tensor, we immediately conclude that 

Ar(&^s + 7) oc -^sin(2mnp), (29) 

precisely in the same way as in Eq. (P6f). The parton expression for T(b — > s + 7) 
scales as 

r (& — > s + 7) oc m 3 Q , rriQ^OQ. (30) 

As a result, our instanton-based model predicts that the duality violating component 
in b — > s + 7 is suppressed as 

AT{b -> 5 + 7)/r (6 -> s + 7) oc — sin(2m Q p) . (31) 

So far it was assumed that the quarks and leptons produced were massless. What 
changes would one decide to take into account finite (nonvanishing) masses of the 
quarks and leptons? 

The exact dependence on the masses of the final quarks or leptons is rather 
sophisticated. Although it is calculable in principle, the corresponding calculations 
are much harder to perform than the simple estimates presented above. Moreover, 
this is hardly necessary. Given a crude nature of the model, which is intended for 
orientation rather than for precision estimates, it seems reasonable to treat u, d and 
s quarks as massless while c as heavy, and the c lines as free propagation. Then the 
only impact of the c quark mass is kinematical, it results in the replacement of the 
total energy itlq by a relevant energy release in the light quarks in the process at 
hand. 

-k -k * 

These are just a few of important applications which are under discussion in 
the current literature. The results are collected in Table 1. Our instanton-based 
model of the duality violation is user-friendly - it is very easy to evaluate the expo- 
nential/oscillating component in other inclusive processes not included in Table 1, 
would such a necessity arise. 

7 Global vs. Local Duality 

Usually by local duality people mean point-by-point comparison of p(s) t heor and 
p{s)ex P , while global duality compares the spectral densities p(s) averaged over some 
ad hoc interval of s, with an ad hoc weight function w(s), 

[■S2 
dsw(s)p(s)theor ~/ ds w(s) p(s) exp . 

S\ J Si 
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H Q -► M 


Hq^X 


H Q ^X 1 


H Q — > H Q ,X 


i 

£3 
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M% 


1 


(am Hq )° 



Table 1: The index of the power fall-off of the oscillating (duality violating) com- 
ponent in various inclusive processes, normalized to the corresponding asymptotic 
(parton model) formulae, in the instanton-based model. The capital X denotes the 
light-quark hadronic states. 



Here I would like to comment on a very common misconception which travels 
from one paper to another. Many authors believe that global duality defined in this 
way has a more solid status than local duality. Some authors go so far as to say 
that while global duality is certainly valid at high energies, this is not necessarily 
the case for local duality. This became a routine statement in the literature. Well, 
routine does not mean correct. 

In fact, both procedures have exactly the same theoretical status. The point- 
by-point comparison, as well as the comparison of p(s)'s (with an ad hoc weight 
function), must be considered as distinct versions of local duality. The distinction 
between the "local" quantities, such as R(e + e~) at a certain value of s and the 
integrals of the type involved, say, in R T is quantitative rather than qualitative. 
Comparison of Eqs. ( |18D and ( p0|) makes this assertion absolutely transparent: in 
both quantities there is a duality violating component, the only distinction is a 
concrete index of the power fall-off (3 vs. 6). 

The genuine global duality applies only to special integrals which can be directly 
expressed through the Euclidean quantities. For instance, if the integration interval 
extends from zero to infinity, and the weight function is exponential, the integral 

dsexp{-s/M 2 }p(s), 



reduces [13] to the Borel transform of the polarization operator II (Q 2 ) in the Eu- 
clidean domain (i.e. at positive Q 2 ). For such quantities, duality cannot be violated, 
by definition. 

There is one more aspect of averaging (smearing) which is not fully understood in 
the literature and requires comment. Let us consider again R T . It can be expressed 
in terms of spectral densities py and pa in the vector and axial-vector channels, 
respectively, 

R r ^ C sf { l - s?) 2 1 1 + 2 w) W') + ^ ■ < 32 > 

(More exactly, the integration over s runs from M 2 , but in the chiral limit we stick 
to, the pion mass vanishes.) The spectral densities py and pA are normalized in 
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such a way that their asymptotic (free-quark) values are 

p v {s) , pa(s) -» N c at s^oo. 

Correspondingly, the asymptotic value of R T is R® = N c . 

In the instanton-based model, the duality violating component in p v< A scales as 
E~ 3 . It is tempting to say then, that the duality violating component in R T can 
be obtained by integrating the duality violating components of pv,A with the weight 
function specified in Eq. fl32l). 



This would lead us nowhere, however. First of all, the integral in Eq. (52[) 
runs all the way down to zero, while Eq. fll8|) is valid at asymptotically large E. 
Even if the lower limit of integration were chosen to be high enough, this would 
not help to find AR T from Eq. (p2|). Indeed, Ap v ^ is an oscillating function of s, 
any smearing inevitably entails cancelations, and the result would depend on subtle 
details of p(s), which are certainly beyond theoretical control. The cancelations are 
seen, in particular, from the fact that asymptotically AR T (M T ) falls off faster than 
Apv,a{s), namely, M~ 6 versus E~ s . There is no way one could predict the change 
of the index from 6 to 3 based solely on the integral ( |32|) . 

At the same time, our instanton-based model does allow one to predict the index 
in R T . To this end one must analyze the appropriate transition amplitude in the 
Euclidean domain performing the analytic continuation to the Minkowski domain 
at the very end. The model captures enough intricate features of QCD to "know" 
the result of the smearing as a whole. Note that AR T (M T ) depends on the highest 
scale in the problem at hand, M T , rather than on any intermediate or low scales one 
encounters in process of integration in Eq. (0). This is a general feature which will 
be valid in any process and will persist in any model based on evaluating singularities 
off the origin. 

The case of R T is quite typical. Similar questions arise in other problems, where a 
similar strategy should be applied. A problem of this type which deserves a special 
mention in view of its practical significance, is that of the inclusive semileptonic 
decays of Hq. The decay rate can be obtained as an integral over appropriate 
kinematic variables over the hadronic structure functions, 



r(H Q ^£u + X) = \V qQ \ 2 -^JdE l dq 2 dq 

{2q 2 Wl + [AEiiqo - E e ) - q 2 ] w 2 + 2q 2 {2E E - q )w 3 } (33) 

where q is the lepton pair momentum, Wi.2,3 are the hadronic structure functions 
which depend on q and q 2 (for their definition see e.g. Ref. p2|). Equation 
(P3|) refers to the rest frame of Hq. Other notations are self-explanatory. The 
exponential/oscillating contribution to the decay rate is presented in Eqs. (|25|) and 
(fffi), where it is assumed that m q = 0. If one decided to get an estimate from Eq. 
(J33|), by integrating the duality violating contributions in the structure functions, 
one would observe parametrically larger violations, exploding at the boundaries of 
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the kinematically allowed domain. These large violations cancel in the total rate 
AT(Hq — >■ h> + X), because of oscillations. 

8 A Resonance-Based Model 

Now we will acquaint ourselves with another approach based on the resonance sat- 
uration of the colorless n-point functions. That this is a distinct dynamical source 
of duality violations follows from the discussion in Sees. 3 and 5. Shortly we will 
confirm this by observing that the functional form of oscillations comes out differ- 
ent compared to that in the instanton-based model. Theoretical analysis is most 
transparent in the limit N c = oo. Later I will allow N c to become finite albeit 
large (Sec. 10). We will first try to abstract general features, and then illustrate 
them in the simplest possible dynamical setting which still exhibits confinement - 
two-dimensional 't Hooft model. f23[ 



Let us first summarize what we expect to get for the polarization operator defined 
in Eq. (fj) in the limit of large N c . In multicolor QCD, N c — > oo, the resonance 
widths vanish. The spectrum of excitations in the given channel is expected be 
(asymptotically) equidistant. A string-like picture of the color confinement naturally 
leads to (approximately) linear Regge trajectories. For each primary trajectory 
there are infinitely many daughter ones. The daughter trajectories are parallel to 
the primary trajectory and are shifted by integers (for a review see e.g. |24j). As a 
result, the excitation spectrum in the given channel takes the form 

Ml = Ml + a 2 n , a 2 = 2/ a' , (34) 

where a' is the slope of the trajectories. Note that the neighboring resonance states 
in the given channel are separated by the interval 2 /at! rather than \joi . This is due 
to the alternating signatures of the daughter trajectories. 

All these properties can be extracted from the Veneziano amplitude found in 
1968 (see Sec. 7.4 in Collins, |24|) which gave rise to the modern string theory. 

Equation (R3j) implies, in turn, that II(g 2 ) can be presented as an infinite sum 

2 N c a 2 °° 1 N c a 2 °° 1 

11(9 } = "12^2 § q 2 ~ Ml = ~Y2^ 2 g 2 - cr% - M 2 ' (35) 

In the problem at hand there are two large quantities: the number of colors and 
the energy. If one fixes s and lets N c grow, one eventually arrives at the comb of 
infinitely narrow 5 functions, as in Fig. 6. On the other hand, if N c is fixed (no 
matter how large it is), with increasing s one eventually finds oneself in the energy 
region where the resonance widths cannot be neglected; in fact, the resonances start 
overlapping, and p(s) gets smoothed. The approximation of the infinitely narrow 
resonances badly fails here, and must be amended. These two limits (N c — > oo with 
s fixed, and N c fixed with s — > oo) are not interchangeable. Later I will include 



the nonvanishing widths (Sec. 10). For the time being I put N c = oo, so that all 
resonance widths vanish. 

The equidistant spectrum only holds if both, the primary and all daughter tra- 
jectories, are exactly linear and parallel. This is not fully realistic. Even putting 
N c = oo does not help. The low-energy parts of the Regge trajectories in QCD 
are not exactly linear, in particular, due to the spontaneous breaking of the chi- 
ral symmetry and the emergence of the massless pions. One can explicitly check 
that the radial excitations are not quite equidistant in the 't Hooft model: they be- 
come equidistant only asymptotically, at n ^> 1 (n is the excitation number). The 
situation in real QCD must be similar. 

This shortcoming of the model affects the condensate expansion at low orders but 
has no impact on deviations from duality at high energies. Since the phenomenon 
under discussion is related to the factorial divergence of the high-order terms, letting 
the low- lying excitations "breathe" does not change the factorial behavior |25||, which 



is in one-to-one correspondence with the spectral formula at large n, where the 
spacings M 2 +1 — M 2 must be constant. Therefore, it is okay to use the linear pattern 
(|35|). For the very same reason the residues of all resonances in Eq. fl35|) are taken to 
be equal. Fluctuations of the residues would show up in the condensate expansion; 
they do not affect the estimate of the duality violations. (For an additional remark 
on the residues see Sec. 11.) 

The infinite sum in Eq. ([33]) reduces to a well-known Euler's ip function, the 
logarithmic derivative of T, 



w ^ z) + ] ' z = — ^ 



U (Q 2 ) = -^ [V»(z) + Const] , : =s V \ ° , (36) 



(see e.g. [ftql ). The irrelevant (subtraction) constant on the right-hand side is infi- 
nite, strictly speaking. The occurrence of the T function reminds us of Veneziano's 
amplitude. 

At positive values of Q 2 an asymptotic representation exists for the if) function, 

^(z) = ln*- — - V^~ 2fc , (37) 

2z j£j 2Jfc V ; 

where B 2 k stand for the Bernoulli numbers, 

a» = (-i)~§P«»>; P8) 

here ( is the Riemann function. (In some textbooks (— l) k+l B2k is called the fe-th 
Bernoulli number and is denoted by Bk-) Equation (|^) defines the asymptotic 
expansion of the polarization operator, 

AT °° H 
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where the coefficients Ck can be expressed through B 2 k (see Eq. (0)) in a relatively 
straightforward manner. The leading logarithmic term exactly coincides with the 
free quark loop of Fig. 1 presented in Eq. (|5|). This explains our choice of the 
resonance residues. The next-to-leading term is 1/Q 2 , followed by higher power 
corrections. Equations Q57|) and (R3) highlight the the factorial divergence of the 



condensate series which I have already mentioned several times. 

One might want to eliminate the 1/Q 2 term from H(Q 2 ) to make the 1/Q 2 



expansion realistic - it is known |14[] that, in QCD with massless quarks, the terms 
of the first order in 1/Q 2 do not appear in U(Q 2 ). The power series starts from 
(G 2 )/Q 4 where (G 2 ) is the gluon condensate. One could eliminate 1/Q 2 , say, by 
fine-tuning the parameter Mq. If Mq = a 2 / 2 = I /a', the 1/Q 2 term cancels. 
Although this might seem desirable, in fact this is hardly worth the bother because 
the model with exactly linear trajectories is too rigid to be realistic anyway. The 
1/Q 4 correction will come out way too large. To accommodate the gluon and mixed 
condensates properly one would need a more flexible model, with more than one 
adjustable parameter. Therefore, I will feel free to simplify the model further by 
putting Mq = and discarding the term with n = in Eq. (pop. 

The spectral density corresponding to the infinite sum of the equidistant reso- 
nances is shown in Fig. 6, 

00 / s \ 2 

Pv ( s ) = PA ( s ) = N c -J2S(--n); a 2 = - . (40) 

n=l w / a 



Let us truncate the power expansion fl39[ ) at some finite order, and examine 
the theoretical prediction for p(s) obtained from the Euclidean side. It does not 
matter in which order we truncate. Any power term (1/Q 2 ) n is invisible in p(s) 
at positive s: analytically continuing to positive q 2 and taking the imaginary part 
one ends up with the S function and its derivatives. The only imaginary part at 
positive s comes from the analytic continuation of InQ 2 . Thus, in the model at 
hand p(s)theor = 1- Comparing this with the comb of the 5 functions in Fig. 6 we 
conclude that the point-by-point duality is maximally violated: between the reso- 
nances the spectral density is grossly overestimated (the "experimental" curve runs 
below the "theoretical" expectation) while at the peaks it is grossly underestimated 
(the "experimental" curve runs above the "theoretical" expectation). What is most 
crucial, the deviations from duality do not die off with energy. The power index 
vanishes! 

All this is pretty trivial. A somewhat less trivial question worth examining is 
the impact of the ad hoc smearing. For instance, let us average the spectral density 
in Fig. 6 with the weight function appropriate to R T , 

I (M 2 ) h(M 2 ) I 3 (M 2 ) 

^-^r- 3 ^r +2 ^T' (41) 

where the moments I n are defined as 

I n (M) = dss n [ Pv (s) + p A (s)} . (42) 

Jo 
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Figure 6: The spectral density in the resonance model. For clarity I gave a tiny 
width to the 5 functions. 

To estimate the oscillating contribution to R T which constitutes duality violation 
that cannot be seen in a truncated OPE we treat M T as a free (large) parameter. 
It will be seen momentarily that for N c = 00 and M T large, yet finite, the duality 
violation in R T scales as 1/M®. In other words, the vanishing power index in p 
translates^ in rj = 6 in R T . Certainly, this distinction is quantitative rather than 
qualitative, but, sure enough, it is quite important from the practical side. 

The sum over resonances in R T is easily calculated analytically: for the spectral 
density of Eq. ([|0) it is 



R T 


= R? 


1E + Ai?° sc 




i?° PE 


= 1- 


a 2 1 
M r 2 + 30 


(a*V 


N c 


\W) ' 



AR°J C 



N r 



a 



2 \ 3 



x(l-x)(l-2x) ( — I + 



*-(!-*)' "30. 
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m • < 43) 



where 



(M 2 \ 
x = fractional part of — j- I , 1 6 [0, 1) 



I presented the result as a sum of two functions of M 2 . The first one, _R° PE , is 
a smooth function expandable in 1/M 2 . The second one, Ai?° sc , oscillates with the 
period er 2 ; its average vanishes, see the plot of AR° SC /R® in Fig. [?|. Although AR° SC 



3 I remind that the power index r\ was defined in Sec. 
M-v sin Mil a 2 . 



3; in the case at hand AR T /R T 
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Figure 7: Oscillations in R T . The plot of AR° SC /R® is presented as a function of 

M 2 r /a 2 . 



is not a pure sine - it contains higher harmonics - the coefficients of the higher 
harmonics are numerically suppressed. With the accuracy of a few percent one can 

write 

Ai?° sc 1 f a 2 \ 3 . (_ Ml\ 

(44) 



li T 




sin 



2vr- 



a 



I pause here to make a comment. The contribution of any particular resonance 
of mass My. to R T , according to Eq. (0), is given by a simple polynomial in 1/Mj? 
(times the step function 8(M 2 — Ml) ). Variations of parameters of a given resonance 
(or resonances) change only the regular terms of the 1/M 2 expansion, but have no 
impact on the oscillatory component. From Eq. (|32f) it is clear that such variations 
change only coefficients of the 1/M 2 , 1/M® and 1/M® terms. OPE for R T must 
exactly reproduce these three expansion coefficients, and so it does. 

In fact, it is not difficult to demonstrate that _R° PE coincides with the OPE 
prediction in the model at hand. The power corrections can be presented as follows: 



* r = Nc+ h_ a b +a h 



(45) 



where the "condensates" J„ are 

/„ ----- / dss 11 [p v (s)+p A (s)-2N c 



(46) 



These integral representations for the "condensates" /„ follow from Eqs. (|32|), ( J42|) 
if one assumes that the spectral densities approach their asymptotic limits faster 
than any power of 1/s. In the model at hand, with the comb-like spectral density, 
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the integral representation (pE6| ) requires regularization. As a regularization one can 
introduce the weight factor exp(— es), taking the limit e — > at the end. With this 
regularization, R° PE from Eq. (j£|) is reproduced. 

The same result for the coefficients of the power terms in R° PE could be obtained 
directly from the expansion of H(Q 2 ) in Eq. fl3"E|). We agreed to put Mq = for 
simplicity. Then the expansion coefficients of II (Q 2 ) are those of the ip function in 
Eq. (|37|) . To find the power terms in R T one may analytically continue the power 
terms in Eq. (|3T| ) to Minkowski, take the imaginary part and convolute with the 
weight function presented in Eq. d32|) . It is easy to see that the only relevant terms 
in Eq. (|37| ) are 1/z and 1/z 4 giving rise to S(s) and S"'(s) in the imaginary part. 
The terms of the higher order in 1/z drop out because the weight function is a 
polynomial of the third order; it contains no s 4 or higher terms. The term 1/z 2 in 
Eq. (^) (it would generate 6'(s)) drops out because the weight function does not 
contain terms linear in s. Then, the expansion for i?° PE is immediately recovered 
provided that one substitutes the appropriate value for the appropriate Bernoulli 
number, £? 2 = —1/30. 

A few words on the numerical aspect. Our consideration is admittedly illustra- 
tive. One should not take too literally the numbers which ensue for many reasons: 
in particular, M 2 t is not much larger than the spacing between the resonances, 
N c = 3 is probably not large enough to warrant the zero width approximation, and 
so on. I would not put too much confidence on particular numbers. I would settle 
on the statement that the above estimate of the oscillation component is valid, say, 
up to a factor of two or so. Taking our formula for AR° SC at its face value and using 
the actual value of the r mass (M 2 /a 2 ~ 1.5) we obtain that AR° SC /R T ~ 3%. It is 
rather rewarding to see that this estimate is in the same ball park as that obtained 
in the instanton-based model. It seems safe to conclude that duality violations in 
the total hadronic r width are expected at a level of a few percent. We do not 
know whether two mechanisms add constructively or destructively. Given this ad- 
ditional uncertainty, it will be no exaggeration to assert that the overall theoretical 
uncertainty 

AR T /R T = 3% . 

The 3% uncertainty in the hadronic r width translates into ~20% uncertainty in 
a s (M T ), which entails the uncertainty of about 6% in a s (Mz). 
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In summary: 

• The resonance structure of the hadronic spectrum associated with the confining 
properties of QCD leads to a distinct exponential/oscillating component invisible in 
the truncated OPE. 

• This component is related to singularities of the appropriate n-point functions 
at infinite separations (in the coordinate space). 
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• The corresponding duality violations are maximal at N c = oo, when the reso- 
nance widths vanish. The power index which was introduced in Sec. 5 is calculable 
in many instances. Generally speaking, the power indices obtained in the instanton- 
based and resonance-based models do not coincide with each other. R T seems to be 
an exceptional case. We do not understand the reasons explaining the coincidence 
of the power indices in R T ; probably, it is accidental. 

• Inclusion of the nonvanishing resonance widths will replace the power suppres- 
sion in the pre-exponent by a weak exponential suppression, see Sec. 10. The change 
of the regimes will probably have little numerical impact in the r decays since ill 2 
is not large enough for the exponential regime to develop in earnest. The factor 
(er 2 /M 2 ) 3 in Eq. (glj) will be replaced by 

/ 2tt5M 2 \ 
6XP {-^^) • 
Both factors are rather close numerically. 

9 Numerical Illustrations in the 't Hooft Model 

The general pattern established in Sec. 8 is nicely illustrated by numerical cal- 
culation in the 't Hooft model which is ideally suited for exhibiting the resonance 
mechanism of duality violations. Indeed, in the 't Hooft model []23| (two-dimensional 
quantum chromodynamics considered in the limit N c — > oo, while a s N c = const) 
the color confinement is automatic because the Coulomb potential in 1 + 1 dimen- 
sions grows linearly with distance. For high excitation numbers the meson spectrum 
grows linearly, M 2 oc n. In two dimensions the gluon field is a nondynamical degree 
of freedom, there are no transverse gluons, and the only remnant of the gluon field 
is the Coulomb (instantaneous) potential. There are no instantons, therefore, the 
mechanism of Sec. 6 does not overshadow the resonance mechanism. Moreover, 
the quark loops are suppressed at N c — > oo - hence, all quark mesons are infinitely 
narrow. The spectrum is calculable, albeit numerically, and so are exclusive decay 
amplitudes, which can be then summed up, one by one. This replaces the experi- 
mental data, to be compared with the inclusive OPE-based calculations. In a sense, 
this is even better than real data in actual QCD which are always incomplete and 
imprecise. In the 't Hooft model one can make gedanken experiments as complete 
and precise as one wishes. 

A thorough analysis of the quark-hadron duality in the semileptonic heavy flavor 
decays in the 't Hooft model was carried out by Lebed and Uraltsev [TIJ. Below I 
will present their result, but at first let me remind relevant aspects of the 't Hooft 
model. 

The dynamical mass scale in the model is set up by the gauge coupling constant 
g (which is dimensionful in D — 2), 

/3 2 ee f(N c - l/N c ) . (47) 
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Thus, (3 plays the role of Aqcd- Being finite at N c — > oo, it provides a natural unit 
of mass. The quarks with masses less than j3 can be called light, while if rriQ 3> (3 we 
are dealing with a heavy quark. To make contact with real QCD the heavy quark 
will be referred to as the b quark, while the corresponding bq bound state as the B 
meson. 

The semileptonic widths b — ► qev are induced by the weak decay Lagrangian 

£ we ak = --^ (q^b) {erfv) , (48) 

where G is a "Fermi constant," and it is assumed that 

m e = m u = . (49) 

Note that in D = 2 the axial current is not independent, so that £ wea k is chosen to 
be pure vector times vector. Moreover, the invariant mass k 2 of the lepton pair is 
always zero ||, provided the lepton masses vanish, Eq. (P3). This fact is specific for 
2D models; it implies that in two dimensions the semileptonic decays b — > qev are 
equivalent to decays into a single massless pseudoscalar particle (b weakly coupled 
to quarks, 

£ weak = — ^qj^e^dvcj). (50) 

V 27T 

For simplicity it will be assumed that the quark q in the transition b — » qev 
is distinct from the spectator light antiquark in the B meson, so that annihilation 
diagrams are absent. The leading (parton-model) result for the semileptonic width 
r(6 — > qev) scales as m\ (to be compared with m\ in real QCD). A few first terms in 
the operator product expansion for V are rather trivial; it is not difficult to obtain 



G 2 m 2 -m 2 m b fi dx 2 
Tope = -: • -rr- / — <p B {x) . (51) 

47r mft M B jo x 

up to terms terms 0(m^ ). Here x is the b quark momentum fraction of B in light- 
cone coordinates, while ^b{x) is the light-cone wave function. The expectation 
value (mb/Ms) (1/x) in Eq. (EH]) can, in turn, be computed in the form of a l/m& 
expansion. 



Equation (pi ) presents a smooth function of mb, a theoretical part of the duality 
relation. It carries no direct information on the opening of new thresholds each time 
Mb crosses the successive values of M n where M n is the mass of the n-th excited 
state in the qq channel (I will denote the meson itself by the same symbol, M n ). 

The "experiential" part is obtained as a sum over individual exclusive decay 
channels of the type B — >■ M n ev. Sure enough, the "experiential" part includes 
the effects of the kinematic thresholds in the most straightforward manner. The 
spectrum and the light-cone wave functions of the excited light mesons are found, 
numerically, from the 't Hooft equation [p3|j . The concrete calculations I would 



like to quote are done []I(| at m q = 0.56/3. The b quark mass in these calculations 

25 



LU 
D. 



O 




rriu 



Figure 8: Deviations from duality in the semileptonic B decays in the 't Hooft 
model, from Ref. |10]. The b quark mass is measured in the units of /3, the light 



quark masses are set at 0.56/?. 

varies from m& = (3 (unrealistically light, no phase space for excitations at all) up 
to mj = 12/3, when the highest kinematically allowed excited qq meson has the 
excitation number 18. (12/3 corresponds, roughly speaking, to m& = 4.5 GeV in 
actual QCD.) The result P| for 



M, 



exp 



1 OPE 

is shown in Fig. 9, as function of the b quark mass. 

Seeds of oscillations inherent to duality violation are clearly seen. (The deviations 
do not average to zero but rather oscillate around a rapidly dissipating contribution 
which can be attributed to discarded higher-order OPE terms.) Numerically, the 
duality violating component turns out to be extremely small - almost certainly, an 
artifact of the 't Hooft model. What I would like to emphasize is not the numerical 
suppression per se, but a general character of the phenomenon. 

10 The Impact of the Resonance Widths 

As was already mentioned, no matter how large N c is, at sufficiently high energies 
one cannot neglect the resonance widths. One should also remember that in the 
real world the number of colors is three - not very large by any count. One can 
expect that the nonvanishing resonance widths, once the resonances start to overlap, 
provide an additional suppression of the duality violating component. 
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Below I will show that, inside the domain of overlapping, the power regime ( |T0D 
is replaced by an exponential one, see Eq. ([0]), even if the power index vanishes 
(as is the case for the spectral density (|3j) presented in Fig. 6). The basic idea is 
that the nonvanishing resonance widths shift the poles away from the physical cut, 
to unphysical sheets, which automatically results in a smoother imaginary part on 
the physical cut, much closer to that obtained from OPE. 

I return to the equidistant resonances considered in Sec. 8 in the leading order 
in 1/N C , with the intention to include the next-to-leading 1/N C effects. For any 
given excitation n, the widths scale as T n /M n ~ N~ l . Below we will analyze an 
infinite sequence of resonances for which we will need to know the dependence of 
this ratio on the excitation number n. For the time being, however, let us focus on 
an "isolated" resonance. 

In the leading order its contribution to the polarization operator II(g 2 ) is 

q 2 

Un (52) 



g 2 -M 2 ' 

here g\ is the residue. If T n ^ 0, in the vicinity of the pole one can use the Breit- 
Wigner formula which replaces Eq. (|52]) by 

' /; ' (53) 



q 2 _ M 2 + lVnMn ■ 

Strictly speaking, g 2 and M 2 in Eqs. (j52|) and (R3|) are different: in the former the 



residue and the mass must be taken in the leading order in 1/N C while in the latter 
they include 0(1/N C ) corrections. This effect is less important than that due to 
r n 's. In many instances (but not always, see below) the 1/N C shifts of g 2 and M 2 
can be neglected. 

Now comes a crucial assertion. At large n the resonance width must scale as 
r n ~ M n /N c , i.e. the n dependence of r„ is the same as that of M n , the square root 
of n. This behavior was predicted long ago |27j on the basis of a simple qualitative 



picture; much later it was confirmed || by numerical studies in the 't Hooft model. 
This formula can be explained as follows. For highly excited states a quasiclassical 
treatment applies. When a meson is created by a local source, it can be considered, 
quasiclassically, as a pair of (almost free) ultrarelativistic quarks; each of them with 
energy M n /2. These quarks are produced at the origin, and then fly back-to-back, 
creating behind them a flux tube of the chromoelectric field. The length of the 
tube L ~ M n /A 2 where A 2 represents the string tension. The meson decay width is 
determined, to order l/N c , by the probability of producing an extra quark- ant iquark 
pair. Since the pair creation can happen anywhere inside the flux tube, one expects 
thatf] T n ~ LA 2 /N C . Taking into account that L ~ M n /A 2 one arrives at 

T n = ^-M n , (54) 



4 Let me note in passing that the 1/N% corrections due to creation of two quark pairs are of 
order L 2 /N% within this picture. Since L ~ M n ~ y/n, the expansion parameter is \fnjN c . 

27 



where B is a dimensionless coefficient of order one. Naively extrapolating Eq. ( |54| ) 
to n = 1 and N c = 3 (well beyond the limits of its applicability) and normalizing by 
the p meson for which Y/M ~ 0.2, one can make an educated guess, 

5-0.5. (55) 

In the 't Hooft model, in which both the meson masses and widths are calculable, 
one can explicitly test the formula T n ~ ^Jn. It is curious that in the 't Hooft mode 
one arrives at a close numerical estimate for B. Figure 9 shows T n versus n. The 
width of the n-th excitation is computed numerically, by summing up all open two- 
meson decay channels, 

r n = £ T{n -> b + c) , M b + M c <M n . 

b,c 

The necessary three-meson constants are found as appropriate overlap integrals 
of the light-cone wave functions, the calculation of T n extends up to the excitation 
number 500. It is seen that the solid curve yielding the y/n dependence fits well 
the numerical data.f] Since M n ~ y/n, the scaling behavior fl5"ip is confirmed. The 
coefficient B defined in Eq. ([54]) is obtained from the coefficient A of Ref. as 
follows: B = A/tt 3 . Since A ~ 15 (see the Erratum), we recover Eq. (|55|) . 

The Breit-Wigner resonance formula now takes the form (in the Euclidean no- 
tation) 



~9n 



V r" - BM l 

}■£ , ft/rz ■ n 



Q z + M< n -i- 



N r . 



(56) 



One must exercise certain care in using the Breit-Wigner expression for the resonance 
contribution in the polarization operator far away from the pole. It must be adjusted 
in such a way that the analytic properties of II(g 2 ) are not spoiled. Namely, II(g 2 ) 
must remain analytic everywhere in the complex q 2 plane, with a cut along the 
positive real semi- axis. No singularities are allowed on the physical sheet, all poles 
must be shifted to unphysical sheets. This is relatively easy to achieve. To this end 
let us replace Eq. flSED by 



~9l 




Q 2 1-— -ln^V +M, 



(57) 



On the upper side of the physical cut near the pole (i.e. at Q 2 in the vicinity of 
— M 2 — ie) the imaginary parts of both expressions (|56| ) and (|57|) coincide, provided 
that g 2 and M 2 are appropriately adjusted (g 2 and M 2 differ from g 2 and M 2 by a 
1/N C correction; in what follows we will omit the tilde, making the adjustment at 

5 Seemingly random fluctuations of T n /y/n around a constant value are also clearly seen. A 
theory of these fluctuations might be a good exercise. Developing such a theory might be fun. 
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Figure 9: The width of the n-th excited meson in the 't Hooft model, in the units 
of 32vr- 2 iV c - 1 /?, from Ref. 0. Note that M n = irf3^/n at n > 1. The parameter (3 
is defined in Eq. 



the very end). Thus, both expressions are equally legitimate for the Breit-Wigner 
description of the resonances. Moreover, Eq. (0), in turn, can be replaced by 



9l 



Q 2 
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a- 
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The distinction between Eqs. (|57| ) and (^) is of the order 0{1/Nj). We do not 
track such terms anyway. This latter formula has the required analytic properties 
- it is nonsingular everywhere in the complex Q 2 plane, except the cut at negative 
real Q 2 . On the physical sheet of Q 2 , the variable 



z = 



'Q 2 



a- 



1-- 



(59) 



never becomes real and negative, so that the pole singularities in Eq. ( |58D are indeed 
shifted to unphysical sheets. 

Summing over the infinite chain of the equidistant resonances, as in Eq. (^), 
we arrive at 



n(Q 2 



N c a 2 



_^ y 

127T 2 1 - B/(tcN c ) ti 



Q 



2 , Q 2 ' 



B 

' wN c 



+ M 2 



a 
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N c 



12tt 2 1 - B/(nN c 



^{z) + - 



(60) 



where z is defined in Eq. ([59D, and 1 — B/(irN c ) in the denominator reflects the 
adjustment of the residues discussed above (this factor can be established by de- 
manding the correct asymptotic behavior at Q 2 — > oo, cf. Eq. (|])). 

How does this compare with the zero-width approximation of Sec. 8? Formally 
both results for the polarization operator look similar; the only difference^ is in the 
definition of the variable z. In the deep Euclidean domain the power expansion 



of U(Q ) ensues from the asymptotic representation (|37|). The k-th. term of the 



expansion, which in the leading order in 1/N C used to be 

(zero- width approx.) 



R / n-2 \ 2k 



now becomes 



2k \Q 2 



B 2k (<r 2 \ 2k A 2lcB, Q 2 




l + 3T^ + 3F + Obi • (61) 



Ik \Q 2 J \ 7iN c 

The impact of the next-to- leading 1/N C terms on the power expansion is quite in- 
significant. In particular, the (B/N c ) InQ 2 correction mimics the logarithmic anoma- 
lous dimension typical of OPE in QCD. 

At the same time, the 1/N C effects radically affect Imll on the physical cut (the 
spectral density) . Instead of the comb of the 5 functions of Fig. 6, we now get at high 
energies a smooth function, with mild oscillations. The plot of the spectral density 
p(s) = (12ir/N c ) Imll (with II determined by Eqs. (BO) and (p5[)) is presented in 
Fig. 10. 

At asymptotically large s the spectral density tends to unity, plus power correc- 
tions corresponding to OPE plus an oscillating component which cannot be obtained 
in the power expansion. One can readily isolate it by exploiting the well-known re- 
flection property of Euler's function, 

1p(z) H = l/}(— Z) — 7TCOt7T2; . (62) 

z 
On the physical cut (in units of a 2 ) 

*- = - sI "*( 1 + ! |) + (w)' (63) 



and the imaginary part of the left-hand side of Eq. (|62|) is essentially given by 



that of — ircotirz. The imaginary part of ip(—z) gives rise to non-oscillating power 

6 The 1/z term which was absent in Sec. 8 is due to the fact that I put Mq = and start the 
summation from n = 1 rather than from n = 0. These simplifications resulting in the occurrence 
of the 1/z term, are irrelevant in the studies of the duality violating component at high energies. 
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Figure 10: The spectral density for the equidistant resonances with finite widths 
specified in Eq. (^4|). 



corrections s~ 2k in the spectral density corresponding to OPE. Alternatively, these 
OPE corrections can be obtained by a direct analytic continuation to Minkowski of 
the power terms fl61"|). 

As a result, the finite-width resonance model leads us to 



p(s) — > 1 + power corr. + 2 exp ( - _ 2 at ) cos ( — — ) , (64) 

where it is assumed that 
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2ixsB 

^-»l but -«1. (65) 

If s is fixed while N c is set to oo we recover unsuppressed oscillations, as in Sec. 
6. In the opposite limit N c fixed and s large we observe an exponential suppression, 
with a weak exponent proportional to 1/N C . 

A key question one can ask in connection with the above analysis is as follows. 
For the given value of N c , what is the boundary energy marking the onset of the 
exponential suppression? The best I can say at the moment is that this boundary 
energy scales as Sq oc N c . A reliable determination of the proportionality coefficient 
is a task for the future. A naive estimate following from Eq. (|6~5D , 

so ~ ?r-£ ~ 2 GeV 2 at N c = 3 , 
ZnJJ 

if correct, would mean that the resonances essentially overlap (and, thus, smear the 
spectral density), starting from the first or second excitation. 
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Figure 11: Experimental data (ALEPH) on the spectral density in the vector channel 
and the duality violation models. 



11 How Does All This Match Experiment? 

In spite of the current extremely high demand on theoretical estimates of duality 
violations, there is surprisingly little effort to elucidate the issue by direct exper- 
imental studies. A high-precision measurement of the ratio R in e + e~ in a broad 
range of energies, from threshold up to, say, s = 10 GeV 2 in a dedicated experiment 
with the proper absolute normalization would give an enormous boost to this issue. 
Alas, such measurements have never been undertaken... The most accurate data on 
the spectral densities were obtained in r decays |p8[ , where they (naturally) extend 
only up to s = 3 GeV . This energy range is way too narrow to be helpful in per- 
fecting the models which are currently in use or in the design of new models. Still, 
confronting the data with current theoretic ideas might give a feeling of whether or 
not we are moving in the right direction. The comparison is presented in Fig. 11. I 
will first explain what is depicted, and then offer several comments. 

The ALEPH experimental data |28|] I show correspond to the spectral density 
in the vector isovector channel. There exist some data above 3 GeV 2 from e + e~ 
annihilation, but the error bars are so large that plotting these points would just 
obscure the picture. Directly measurable in the hadronic r decays is the sum of 
the vector and axial spectral densities. To obtain the spectral densities separately 
one has to sort out all decays by assigning specific quantum numbers to each given 
final hadronic state. In the majority of cases such an assignment is unambiguous. 
For instance, two pions (whose contribution is the largest) can be produced only 
by the vector current. Some processes, however, can occur in both channels, for 
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instance, the KKtt production. Using certain theoretical arguments it was decided 
28] that around 3/4 of all KKtt yield must be ascribed to the vector channel. Other 



theoretical arguments El , which seem more convincing, tell that virtually all KKtt 



production takes place in the axial channel. Therefore, in fitting the data, I subtract 
the KKtt yield from the data points presented in Fig. 11. I hasten to add that the 
subtraction affects only the energy range s > 2 GeV 2 , i.e. the second peak in Fig. 
11, and even in this energy range the effect is small, ^5%. The subtraction is not 
essential for a general picture I draw here. 

The solid curve at s < 2.7 GeV 2 is the best fit of the data points thus obtained 
by the sum of two Breit-Wigner peaks (the first one, the p meson, is a modified 
Breit-Wigner taking into account threshold effects important for the p meson). For 
my purposes the two-resonance fit presents an excellent approximation to the ex- 
perimental spectral density at s < 2.7 GeV 2 . 

The tails of the curves at 2.7 < s < 7 GeV 2 represent the resonance-based and 
instanton-based models (the solid and dashed curves, respectively) which I discussed 
in Sees. 10 and 6. The solid curve is a modulation of p(s) p0 rt by the factor 

1 + 1.22E' 3 sin(2pE-S), 

p = 3GeV~ 1 , 5 = 1.32, E in GeV, 
cf. Eq. flilf) . The dashed curve is a modulation of p(s) per t byQ 

/ 2TrsB\ (2tTS 

1 - L24exp (-^ si,1 (^- 308 

a 2 = 2 GeV 2 , 5 = 0.5, iV c = 3 , 

cf. Eq. (H). 

Finally, the thick solid curve in Fig. 11 displays p(s) pcrt calculated through order 
a 3 s , with A® = 200 MeV. 

IVlo 

It is clearly seen that, qualitatively, both models match the data well. In fact, in 
actuality I would expect an oscillating component which is a combination of the solid 
and dashed curves, since the mechanisms of the duality violations they represent are 
complementary. The instanton mechanism leads to a visibly slower (power) fall off, 
with rearer oscillations. It is also clear that experimental measurements at the 
percent level of accuracy will most certainly provide us with the material needed for 
construction of a reliable and well calibrated model of the duality violations. 



7 As we know from Sec. 10, the equidistant resonances with equal residues result in oscillations 
of this type superimposed on/?= 1. By a trivial adjustment of the residues one can readily achieve 
that the oscillations (p3) are superimposed on p(s) pe it- 



33 



12 "Exclusive" Duality 

This section is somewhat perpendicular to the main theme of the review, and can 
be safely omitted. While previously my topic was estimating duality violations, 
now I will discuss a special "exclusive" mode of implementation of the quark-hadron 
duality. 

In all problems considered above duality is applicable at high energies (momen- 
tum transfers), where the cross sections (decay probabilities) are saturated by a 
large number of exclusive channels. Moving in the opposite direction, towards lower 
energies (momentum transfers), we decrease the number of the open channels and, 
typically, worsen the accuracy of the quark-hadron duality. There is a stereotype in 
the public eye that "duality works well in inclusive processes after summing up over 
a large number of exclusive channels." Although this is often true, sometimes the 
stereotype turns out to be wrong. In this section, concluding my review, I give two 
examples when duality (i.e. "equality-with-controlled-accuracy" between the quark 
and hadronic processes) takes place even though the hadronic process is saturated 
by a single exclusive channel. Sure enough, such situations are rather exceptional 
and are explained by certain custodian symmetries. 



The first example of this type was discovered |30| in the mid- 1980 's in the heavy 



quark transitions at zero recoil. A prototype process has the form 

Q -> Q' + tv (66) 

or, at the hadron level, 

H Q ^H Q ,+£v. (67) 

Here Q and Q' are heavy quarks of distinct flavors, Hq stands for a Q-containing 
heavy hadron (practically, the lowest-lying state, say Hq = B or Hq = A b ), while 
iv is the lepton pair. Needless to say that the sum over all possible Hq/'s in the 
final state is implied since we consider the inclusive reaction. As we will see shortly, 
at zero recoil, only one state will survive in this sum. 

At the point of zero recoil, the total spatial momentum of the lepton pair van- 
ishes, k — 0. Then the time component of the lepton pair momentum in the physical 
decay (|67|) must be equal to k = Mh q — Mu . The value of ko is maximal if Hqi 
is the ground state. It is convenient to introduce the notation 



AM=(M„ g -M H J dM> (« 



Q / ground state 

and 

e = AM - k . (69) 

Then on the physical cut e is real and positive. 

Dynamical information is encoded in the differential probabilities which go under 
the name of the hadronic structure functions Wi (they were already mentioned in the 
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Figure 12: Inclusive semileptonic decay of one heavy flavor into another at zero 
recoil, (a) Free quark transition operator; (b) hadronic saturation. The closed circle 
denotes the insertion of the axial or vector current. The four-momentum of the 
lepton pair is denoted by k^. 

very end of Sec. 7). One starts from the transition operator describing the forward 
scattering of Hq to Hq via intermediate states Hq>, 

% v = iJd 4 xe- lkx T{Q(x)^ 5 Q'(x) , Q'(0)^ 5 Q(0)} . (70) 

For definiteness I limit myself here to the axial- vector current, Q^^zQ' ■ (The vector 
current can be treated in a similar fashion.) The hadronic amplitude is obtained by 
averaging T^ u over H Q , 

Ku = ^(H Q \f, u \H Q ) . (71) 

From kinematics we infer a general decomposition 

/V = -hig,u, + h 2 v ll v v - ihaeuwapVakp + h 4 k fl k v + h(k^v u + k v v^ , (72) 

where v^ is the four-velocity of Hq (in the rest frame t> M = {1,0,0,0}.) At zero 
recoil h^s depend only on fc . The physically observable quantities are the structure 
functions 

wt = 2 Im fn , < e < AM . 

The relation between h^s and w^s is the same as that between H(Q 2 ) and p(s) 
in e + e~ annihilation. Much in the same way as in the latter problem, one can 
calculate hi far away from the physical cut, then perform analytic continuation onto 
real positive e, and then take the imaginary part, term by term in OPE. A new 
element in the heavy flavor decay is that the expansion runs now in two parameters, 
A/e and A/rriQ, where tuq is the heavy quark mass. It is assumed that K/uiq^qi <^ 1. 
Shortly I will send the quark masses itlqqi to infinity keeping the difference vriQ—mQi 
fixed. 

At the quark level h^s are determined by the diagram displayed in Fig. 12(a). 
By analogy with e + e~ annihilation one might think that, to make OPE meaningful, 
it is necessary to require A/e <^i 1. Surprisingly, at zero recoil this hasty conclusion 
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is wrong. In fact, it is not difficult to prove (see e.g. f5If ) that the structure of the 
power expansion in the case at hand is as follows 

-h = r( 1 + Ar + -) 

™Q - m Q > - fc V m Q J 

A 




[m Q - m QI - 
+ 7 ^ . v , l£r + -l+-- (73) 

[m Q - m QI - k Q ) 
If one takes the limit tuq — > oo first, keeping finite the quark mass difference, 

Am = rriQ — tuqi = fixed , (74) 

then the expansion for hi becomes absolutely trivial, namely 

- hi = (75) 

mQ - m Q > - k 

exactly, for all values of tuq — vhqi — ko, large or small. Thus, at zero recoil, the 
free quark amplitude of Fig. 12(a) is exact, a manifestation of the heavy quark 
symmetry. 

Let us now examine hi from the hadronic side (Fig. 12(b)). The very same 
heavy quark symmetry tells us that the only state surviving in the sum over Hqi at 
k — 0, rag' = oo is the ground state (Q'q) where q stands for the light spectator 
quark. Moreover, 

F 2 

- (^l)hadr = , (76) 

e 
where F is the Hq — > Hqi form factor at zero recoil. Since AM = Am, both 
expressions (|75|) and (76) perfectly match provided the transition form factor is 



unity at zero recoil |5D[. The quark-hadron duality is perfect: the quark pole of Fig. 
12(a) is exactly equal to the hadronic pole presented in Fig. 12(b). Correspondingly, 
the same equality takes place for the structure functions Wi(ko). 

The second example of "exclusive" duality I would like to mention is somewhat 
more exotic. It was found || as a byproduct in the studies of the Pauli interference 
in the heavy flavor decays in the 't Hooft model. The setting has already been 
presented in Sec. 9, which is devoted to the semileptonic decays. Now I pass to 
nonleptonic decays. Correspondingly, the weak Lagrangian in Eq. ( |4*8| ) is replaced 
by 

£weak = —7= {oi (cj„b) (d-fu) + a 2 (d^b) {crfu)} + H.o, (77) 

where I use the same notation as in Sec. 9; two distinct four-fermion operators in 
Eq. (|77|) , with (dimensionless) coefficients ai and a 2 , represent two possible patterns 
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Figure 13: The quark graph determining the effect of the Pauli interference in the 
nonleptonic inclusive B~ decays at leading order in OPE. 

of the color flow, direct and twisted (in actual QCD the latter emerges due to hard 
gluon exchanges). At the quark level the Pauli interference in the B~ decays is 
described by the graph of Fig. 13, which determines the coefficient in front of 
the four-fermion operator (lyy^u) (uj^^b) in OPE. To calculate the effect of the 
Pauli interference in the quark language one has to calculate the diagram, take the 
imaginary part and average the four-fermion operator over B~, 



2M F 



(B~ 



(h,l 5 u) (uYl 5 b) | B-) = l -f B M B . (7 



If nib ^> (3, the result thus obtained has the same status as the parton-model formula, 
say, for R(e + e~). 

Following Ref. 9, I will set all quark masses except that of b to zero, m c = m u = 
Tad — 0. Then the Pauli interference of Fig. 13 takes especially simple form 

Arpi(B-) = -^(a ia2 )G 2 f B M B . (79) 

This effect is a l/m& correction to the total inclusive B~ width.0 At the same time, 
AFpi is of the same order in 1/N C as r tot (I? _ ). This is important. 

As far as the hadronic saturation is concerned, the following statements are valid 
to the leading order in 1/N C : 

(i) The intermediate states saturating AFpi are, by necessity, two-meson states. 
Examples are shown in Fig. 14. 

(ii) The graph of Fig. 14(a) is due to the first operator in Eq. (|77| ) and is 
proportional to ai, while that of Fig. 14(b) is due to the second operator and is 
proportional to a 2 - 

(iii) Factorization applies. 

Generally speaking, both mesons in the intermediate state, ir~ and D°, need 
not be ground states. Radial excitations would be acceptable, too. However, if 
m c — m u — m d — 0, the only hadronic states produced by the currents d'j^u and 



8 In actual QCD the Pauli interference would be a 1/m^ correction; 1/m;, versus 1/mjj reflects 
the fact that the canonic dimension of the current q^^q is one in D = 2 and three in D = 4. 
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Figure 14: The product of these two amplitudes determines Ar PI to the leading 
order in 1/N C . 

u'y^c are the ground state pseudoscalar mesons with zero mass. This is a special 
feature of the two-dimensional theory, with no parallel in four dimensions. The fact 
that (du) meson in Fig. 14(a) is produced by the dj^u current forces it to be a 
massless tc~ , while, by the same token, the (uc) meson in Fig. 14(b) is a massless 
D°. 

As a result, the only intermediate hadronic state surviving in the Pauli interfer- 
ence is 

B~ —> massless tt~ + massless D° — > B~ (80) 

This single exclusive channel must saturate (|7^) . 

So it does! This comes out rather trivially || because of the following properties 
of the pion coupling and the B —>■ D(tt) transition form factors at k 2 = 




k^-±-(n-\d lu b\B-) = -k z f cj> B {x) dx = -k z f B 

2Mb Jo V ly r . 




(81) 



Here (j)(x) is the light-cone wave function; for the pion <p(x) = 1. Note the occurrence 
of Jb in the transition form factor. 

Using Eq. (J3l"l) it is easy to verify that the massless quark loop of Fig. 14(a) 
perfectly matches the massless meson loop, the only contribution to AFpi surviving 
at the hadronic level. The origin of the factor f B in Eq. (|79f ) in the hadronic 
calculation is totally different, though; it comes from the square of the transition 
form factor, the matching seems magic, yet it could not have happened otherwise. 

It is worth stressing that the magic simplifications which led to the "exclusive" 
mode of the duality implementation in the problem at hand are explained by mirac- 
ulous features of two dimensions. 
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13 Conclusions 

In this review I identified general mechanisms causing deviations from duality and 
derived scaling laws governing the damping of the duality violating component in a 
variety of inclusive processes as a function of energy (momentum transfer). The esti- 
mates of the absolute normalization one can perform at the moment are less certain. 
This explains why I was so cautious in my discussion of the numerical situation in 
R T , and completely avoided this issue in other inclusive processes. Getting a better 
idea of the absolute normalization is absolutely necessary for all practical applica- 
tions of the theoretical constructions presented here. Precision measurements of 
p(s) in a wide energy range would be of enormous help in this question and would, 
probably lead to a speedy solution. 

As I have already noted, duality violations parametrize our ignorance. Were an 
analytic solution of QCD found, the contents of this review would become instantly 
obsolete. Then the exact asymptotic behavior of inclusive cross sections would be 
known, and the very concept of duality violations would become irrelevant. 
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